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REFERENCES TO RELATED APPLICATIONS 
[0001] This application is a continuation-in-part of United States Patent Application 
Ser. No. 09/949,966 filed on September 10, 2001 that is a continuation-in-part of 
United States Patent Application Ser. No. 09/542,307, now United States Patent 
******** ? Q f are incorporated herein by reference. 

BACKGROUND OF THE INVENTION 
Field of the Invention 

[0001] The invention relates to a method for modeling stresses in the vicinity of a 
borehole and predicting the failure of liners or screens. 

Background of the invention 

[0002] Sand control screens are utilized for various purposes in subterranean wells. 
The name derives from their early use in preventing the production of sand along with 
fluids from formations. A sand control screen is typically suspended from production 
tubing extending to the earth's surface and positioned in a wellbore opposite a 
productive formation. The wellbore in an annular area between the screen and the 
casing may be filled with a relatively large grain sand ("gravel"). This gravel 
prevents the fine sand from packing around the production tubing and screen and the 
screen prevents the large grain sand, gravel, from entering the production tubing. In 
this way, the sand control screen may exclude the produced sand while permitting the 
valuable fluids to enter the tubing for transport to the earth's surface. Perforated or 
slotted liners or expandable liners may also be used instead of a screen as a sand 
exclusion device. Hereafter, for the purposes of this application and the description 
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of the invention, the term liner is used to refer to all such types of sand exclusion 
devices. 

[0003] There are numerous prior art devices with different configurations of the liner 
5 and the components of the gravel pack and numerous methods for deployment of the 

devices. Examples of such devices are found in US Patent 5,829,522 to Ross et aL, 
and US Patent 6,053,250 to Echols. Regardless of the device used and the method of 
deployment, the line and gravel pack must be designed to withstand the stresses in 
the subsurface formation: failure of the screen or a breakdown of the gravel pack can 
1 0 lead to sand production from the formation and plugging of the borehole. The 

selection and design of the liner depends on analysis of formation strength, strength 
distribution, permeability, permeability distribution, shale content, fines migration, 
and grain size and distribution. Of these factors, permeability, permeability 
distribution, grain size, and grain distribution can be measured with reasonable 
15 accuracy; however, most oil companies still have difficulty conducting formation 

strength and strength-distribution analyses. The two primary reasons for poor 
analysis are that mechanical logs available from service companies are not reliable or 
must be calibrated with other methods and that reasonably reliable numerical models 
for strength analysis are owned exclusively by several companies. 

20 

[0004] Another factor complicating the design of liners is that when fluid is produced 
from a reservoir, a reduction in pore pressure occurs. The pressure reduction is 
greater near the wellbore and increases with production time and rate. The reduction 
in pore pressure causes compaction of the formation containing the reservoir fluid, 
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which imposes radial and axial loads on the well. Wellbore loads resulting from 
reservoir compaction are seldom considered in the design of casings, liners, and 
gravel-pack screens, yet they can be significant. Radial and axial pressures on 
wellbore tubulars from reservoir compaction are illustrated in Fig. 1 for a vertical 
well. For a deviated well, the wellbore is exposed to a radial component of the 
vertical overburden load, increasing external pressures on tubulars. Determining 
reservoir compaction loads on wellbore tubulars is not a simple task. Field 
measurement of reservoir compaction loads is difficult because of the time required 
for these loads to develop and the difficulty in measuring them. Simple analytic 
techniques for calculating reservoir compaction do not account for all the important 
variables affecting well loads. A computer model, however, can incorporate the 
many important variables to obtain realistic predictions of these loads. 

[0005] Morita provides a comprehensive discussion of different sand-prediction 
models. As discussed therein, three types of field models are commonly used for 
predicting the onset of sand production. Core-based models are numerical models 
with a set of stress-strain curves as input to calculate the cavity stability and strength. 
Sonic based models are direct approximations from standard mechanical logs that 
calculate cavity stability using a linear stress/strain model. Regional statistical sand- 
production models involve back-calculation methods where a base solution 
incorporates a correction factor determined by field sand-production data. 

[0006] Wooley & Prachner provide a methodology for the analysis of compaction 
loads on casings and liners of vertical and/or deviated boreholes resulting from 
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reservoir depletion. The methodology uses a finite element simulation and is based 
on the separation of formation stress into a component of pore pressure and matrix 
stress. First, overall subsidence is computed with a large-scale model. The solution 
applies undisturbed boundary conditions far from the compacting reservoir and 
computes deformations and stresses near the well, in the reservoir, and in surrounding 
formations. Second, a near-well analysis is used to compute loads on the well. 
Effects of formations beyond the near- well region are transmitted to the solution by 
means of boundary displacements, which are defined from deformations computed 
with the large-scale model. 

[0007] As discussed by Hamidet al, liner collapse occurs from application of 
differential pressures (across the screens) that exceed the collapse strength of the 
screen jacket. The resulting failure may not always be critical. In general, screen 
jacket collapse may lead to subsequent long term erosive failure. 

[0008] Modeling of liner failure is typically carried out using a finite element analysis 
(FEA). Guinot et al. (US Patent 6,283,214) use a FEA to show that a particular shape 
and orientation of the perforations of a liner minimizes this destabilization, hence also 
minimizes sand production. In particular, and in the specific case of a vertical 
wellbore, for instance, elliptically shaped perforations, having the major axis aligned 
in the direction of maximum principal in situ, or compressive stress, improve the 
stability of the formation in the region near the wellbore, hence minimizing sand 
intrusion. Particularly preferred embodiments of this aspect of the Invention are 
perforations with an aspect ratio of about 5:1, and having their principal axis 
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substantially aligned with the direction of maximum compressive stress. 

[0009] Prior art methods of modeling the failure of liners predict the load on the liners 
using a combination of experimental data and FEA. The experimental data are used 
5 to define parameters of a model in which the parameters characterizing the liner and 

the formation are relatively uniform. There may be a variation in the stresses applied 
to the liner but the material comprising the liner and the formation is relatively 
homogenous. Failure prediction in models like this is based on average properties of 
the formation surrounding the liner as well as the liner itself. In reality, such 

10 homogeneity rarely exists in the formation around the well and in the liners: on a 

local scale, there may be statistical variations in the strength including anisotropy 
effects. In practice, failure usually occurs at the weakest point resulting in 
asymmetric loading creating localized deformation that can cause early failure of the 
liner. Accounting for statistical distribution of properties and the development of 

1 5 local failure areas are impractical to consider using the commonly used finite element 

analysis techniques. There is a need for a method of analysis and prediction of failure 
of liners that takes into account statistical variations in the properties of the vicinity of 
the borehole and liner. Such an invention should also be computationally fast. In 
addition, it is preferable that the invention should be user friendly in that specification 

20 of the material properties and loading be easily input and that the invention be able to 

provide graphical displays of the deformation and fracture process. The present 
invention satisfies this need. 
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SUMMARY OF THE INVENTION 
[0010] The present invention is computer implemented method of modeling failure 
of a borehole in a subsurface formation. A mode of definition of a subsurface model 
characterizing the borehole and its environs is selected that is one of (i) an aerial 
mode wherein the model comprises a plurality of nodes in a plane orthogonal to a 
longitudinal axis of the borehole, and (ii) a 3-D mode. A subsurface model is defined 
that includes the borehole and at least one additional region selected from (i) a liner 
in the borehole, (ii) a casing in the borehole, and (iii) at least one earth formation such 
as a near region and a far region. Each of said plurality of regions comprises a 
plurality of nodes interconnected by linkages that may be springs, rods or beams. 
Material properties associated with the regions are defined, the material properties 
having a statistical variation. An initial deformation pattern of the model is defined 
and using a dynamic range relaxation algorithm (DRRA) a force equilibrium solution 
for the model and the initial deformation pattern is determined that includes 
fracturing that simulates failure of the material surrounding the borehole. In an 
alternated embodiment of the invention, instead of an initial deformation pattern, 
stresses are applied on the boundaries of the model. 

[0011] Stresses applied in the model may include stresses produced by reservoir 
depletion and an associated decrease in formation fluid pressure. This makes it 
possible to predict lining of casing failure in a producing reservoir. 

[0012] In an alternate embodiment of the invention, large scale geologic simulation 
using a prior art DRRA method is used for determining stresses in subsurface 
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formations. Using these stresses, it is possible to use the method of the present 
invention to predict possible failure of a wellbore during the process of drilling along 
a predetermined trajectory in the subsurface. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

[0013] The file of this patent contains at least one drawing executed in color: Copies 
of this patent with color drawing(s) will be provided by the Patent and Trademark 
Office upon request and payment of the necessary fee. For detailed understanding of 
the present invention, reference should be made to the following detailed description 

1 0 of the preferred embodiment, taken in conjunction with the accompanying drawings, 

in which like elements have been given like numerals and wherein: 
Figure 1 (Prior Art) illustrates the phenomenon of reservoir compaction due to 
pressure reduction. 

Figure 2 (Prior Art) is a flow chart illustrating the major steps of the use of a 
1 5 Dynamic Range Relaxation Algorithm for modeling of fractures. 

Figure 3 (Prior Art) show the triangular nodal configuration for an aerial model. 
Figure 4 illustrates the geometry of an illustrative model with four zones. 
Figure 5 shows the nodes for the model of Figure 4. 

Figure 6 shows the nodal configuration of the model of Figure 4 after application of 
20 stress leading to failure of the borehole. 

Figure 7 is a post-failure view of the model emphasizing the fracturing in the model. 
Figure 8 shows post failure views of the model for two different values of the near 
rock compressibility. 

CON- 1021 8 



DETAILED DESCRIPTION OF THE INVENTION 



[0014] The present invention uses a Dynamic Range Relaxation Algorithm (DRRA) 
for the modeling of borehole failure. United States Patent Application Ser. Nos . 
09/542,307 (the '307 application) filed on April 4, 2000, now United States Patent 
**********^ and Application Ser. No. 09/949,966 (the '966 application) filed on 
September 10, 2001 disclose a method of using a DRRA for the modeling of 
deformation and fracturing of earth formations on a geologic scale. The present 
invention uses many of the concepts from the '307 and the '966 application. 

[0015] Turning now to Fig, 2 (Prior art) , a flow chart of the major steps of using a 
DRRA are shown. The first step in the invention is to select a mode of definition of 
the subsurface 101. This step defines the boundaries of the model and the nodal 
configuration therein. The mode of definition may be aerial, cross-sectional or 3-D. 
Within the model, a plurality of interconnected nodes that characterize the geometry 
of the model are defined. In a preferred embodiment of the invention, the nodal 
pattern is a regular triangular lattice, although other patterns, such as a random lattice, 
may also be used. The user may also specify the number of nodes in and the aspect 
ratio of the model. 

[0016] Within the framework of the nodal geometry defined at 101, the material 
properties of model are input 103. The nodes are interconnected by bonds such as 
springs, beams, or rods having elastic properties and breaking strengths related to the 
physical model. In a preferred embodiment of the invention, the springs are linear 
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elastic springs. In an alternate embodiment of the invention, any user-specified 
stress-strain curve may be used. The user may also specify independently a repulsion 
between the nodes. In the beam model, the force is based upon linear beam equations 
of standard elastic theory. In the rod model, there is an angular force determined by 
5 the angle between links between nodes. The purpose of these forces (spring, beam, or 

rod) is to stabilize the matrix involved in the solution of the deformation process. In a 
preferred embodiment of the present invention, the conditioning of the model 107 as 
described in the '307 and ' 966 applications is not used. 



1 0 [0017] In one embodiment of the invention, the forces act between adjacent nodes 

(nearest neighbors). In an alternative embodiment of the invention, the forces act 
between additional nodes that are farther away, i.e., between next nearest neighbors or 
even further neighbors. Again, the addition offerees between next nearest neighbors 
is to stabilize the solution matrix. In a preferred embodiment of the invention, the 

1 5 next nearest neighbor forces are used in conjunction with the spring model, though 

the next nearest neighbor forces could also be used with the beam model or the rod 
model. 



[0018]Once the model has been defined, the deformations are applied to the model 
20 109. The result of deformation is to produce a deformed model with faulting and 

fracturing therein. This determination of the deformation process is carried out using 
a Dynamic Range Relaxation model 115 that is discussed in the 4 3 07 and ' 966 
applications. The resulting deformed structure 117 from the model simulates the 
fracturing processes that occur in the real world and may be analyzed by those versed 
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in the art to determine a likely mode of failure of the borehole. 

[0019] The '307 and '966 applications also discuss use of the Anticipate 111 and 
model adjustment 113. These are intended for getting an approximate look at the 
5 results of the deformation and are primarily intended for the large complicated models 

encountered in geologic modeling. These steps are not used in a preferred 
embodiment of the present invention that is aimed at local modeling of borehole 
failure. They may, however, be used when simulating the full length of a borehole. 

10 [0020] Another aspect of the '307 and '966 inventions is a conditioning process. The 

geologic interpretation process starts with the final position of the faults, and, in 
particular, the large-scale faults 105. These could be observations of surface faulting, 
as well as faults interpreted from the wellbore. In geologic modeling, it is essential 
that the deformed structure 117 include in it at least these large scale faults. Usually, 

1 5 the initial model is obtained from basic principles that would be familiar to a 

structural geologist. The observed fault structure is "undeformed" by reversing the 
process that produced the faulting in the first place. This is entirely a kinematic 
procedure that repositions the various fault blocks consistent with the laws of gravity 
and conservation of mass to the initial position they must have been in before the 

20 deformation started. This process is sometimes referred to as palinspastic 

reconstruction. Application of deformation to the FramView™ model used in the 
'307 and '966 applications, usually results in the observed large-scale deformation. 
However, in some cases this may not be enough. In such cases, the conditioning of 
the model is carried out. In the conditioning process, the model is weakened at the 
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reconstructed fault positions so that after deformation of the model, fracturing and 
faulting are more likely at these positions. This too is not usually necessary for local 
borehole failure modeling. 

[0021] The '377 and '966 applications also provide a discussion of the graphical user 
interface used in geologic modeling. With slight modifications, the graphical user 
interface may also be used in the present invention. 

[0022]Turning now to Fig. 3, an aerial model is depicted. In the picture above, two 
sets of nodes are shown. The lower plane is a set of nodes 201a, 201b, 201c. . . 201n 
called the "substrate". The upper plane consists of the nodes 211a, 211b, 211c . . . 
21 If that make up the material itself. In one embodiment of the invention, called the 
nearest neighbor model, a line is drawn connecting each of the nodes in the material 
with its neighboring nodes. These lines represent bonds such as springs, beams, rods 
or attractive forces between the nodes. In the example here two sets of springs are 
shown. One set is in the plane of the nodes; the other extends down to one of the 
nodes in the substrate. In a way, this picture is incorrect, because the actual model is 
really only two-dimensional. The entire model lies in a plane, so that the substrate 
nodes are not really below the material nodes, but instead lie in the same plane with 
the material nodes. In an alternative embodiment of the invention, forces are 
introduced not only between nearest neighbors but also between next nearest 
neighbors. 

[0023] To deform the material, the substrate nodes are moved according to 
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deformation arrows defined in the model interface. In the aerial mode, deformation 
arrows (not shown) car be applied at any point in the substrate. These are applied in 
a point and click manner. A suitable interpolation is done between the points at 
which the deformation is applied. As the substrate nodes are moved, the material 

5 nodes are "pulled along" and are separated from each other. As part of the model 

definition, a statistical distribution of breaking strengths is defined in the interface. In 
a spring model, the distribution of breaking strengths has an equivalent distribution of 
breaking lengths. Hereafter, the two terms may be used interchangeably. From this 
distribution of breaking lengths, a breaking distance has been assigned to each spring. 

1 0 If the length of any spring exceeds its breaking distance, that spring is broken and is 

never re-attached. The distribution of breaking lengths is assigned by giving a mean 
and the standard deviation of the breaking distances. A random number seed is used 
to choose one breaking length from the distribution for each spring on the spring 
network. As soon as the length of a spring connecting two nodes exceeds its breaking 

1 5 length, the spring breaks. By changing the random number seed, a number of 

realizations can be made from the same breaking distribution. These features cannot 
be implemented in the usual finite element analysis approach without the user 
recreating the finite element mesh as each fracture develops. This makes the finite 
element approach practically impossible when large numbers of fractures 

20 progressively develop. 

[0024] In the 3-D mode, the nodes are preferably in a triangular configuration with an 
overall cylindrical geometry to simulate a borehole. The deformation may be applied 
to the sides of the volume of the material region. 
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[00251 The FramView™ simulation requires that at each time step the sum of forces 
on each node to be zero. The forces described in the "307 and '966 applications 
consisted of the spring forces (or beam or rod) and gravity. In the present invention, 
another force can be added to directly to each node. In one embodiment of the 
5 invention, this is the difference in fluid pressure between one side of the node and the 

other. This fluid force gradient can be calculated by solving, in addition to the 
discrete element simulation (FramView™), a finite difference (or finite element) 
simulation which tracks fluid flow. These flow simulations go by the name of 
"reservoir flow simulators" and are well known to those versed in the art. These flow 
; A 1 0 simulators calculate the fluid pressure at any point in space in response to fluid 

; 3 production or injection at well bore locations. Knowing the fluid pressure at each 

i H point also means that the pressure gradient is known at each point. In an optional 

| ~; embodiment of the invention, this pressure gradient at each node location can be used 

; A to predict an additional fluid force on each node. With this modification, in addition 

15 to having the fluid simulator influence the node motions, the node motions in turn can 

i 
I 

i J be used to influence the fluid flow simulations. One of the input values in a fluid 

flow simulator is the permeability (or transmissibility) between each point modeled in 
the finite difference or finite element grid. The change of any node positions can be 
used to influence the local permeability in that region of space. For example, if a 
20 bond is broken due to node motions, permeability in the flow simulator can be 

increased to simulate increased ease of flow due to a fracture a that location. 
Similarly, bonds pushed tightly together can cause a decrease in permeability if there 
are no broken bonds. Lastly, if desired, the change in local density of nodes could 
also change the porosity in the flow simulations (which is also an input parameter 
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needed for these simulations.) In prior art flow simulators, permeability and porosity 
are not varied in time, but are assumed constant. In an optional embodiment of the 
present invention, permeability and porosity are changed as nodes moved in the 
discrete element simulation. The invention described in the '307 and the 6 966 
5 applications has been modified in the present invention so that instead of applying 

deformation to the substrate, stresses may be applied to the boundary nodes as 
described earlier in this paragraph. With the ability to add stresses directly to nodes, 
the FramView™ software has additional flexibility in the simulation of deformation 
processes. An example is given next. 

10 

[0026] In an optional embodiment of the invention, the initial geometry may be 
defined by using a FramView™ simulation of large scale geologic deformation such 
as described in the '306 and '977 application. A geologic-scale simulation is run to 
produce a deformation that now appears in the subsurface. For example, one side of a 

1 5 material region is extended until a graben is produced that looks like what is seen in a 

seismic image. This may be done in 2D or 3D. The "final 55 state of this FramView™ 
geologic simulation would have not only the positions of the nodes corresponding to 
material in the subsurface, it would also contain the forces or stresses between the 
nodes. In one application, a well trajectory is defined through the material region 

20 which would correspond to a planned well track. The stresses from the geologic 

simulation are then included in the present invention as preexisting stresses and used 
to identify points along the well track that are likely to result in failure. In an 
alternate embodiment of the invention, nodes along the well track are removed to 
simulate the drilling of the well. Removing these nodes would cause the remaining 
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nodes to move in response to the changing stress which would occur because nodes 
are now "gone". If the nodes surrounding the wellbore tended to move into the space 
occupied by the missing nodes, this would indicate that in these regions the well is 
more likely to fail at these locations rather than other locations along the well track. 
This simulation of wellbore failure does not require the presence of a casing or liner 
in the borehole and the model may just comprise the borehole and the formation. 

[0027] Turning now to Fig. 4, an exemplary model of the vicinity of a borehole is 
shown. Four regions are indicated by the wellbore 301, a liner 303, the near rock 
earth formation 305 and the far rock earth formation 307. The four regions are for 
illustrative purposes only and the method of the present invention may be used for 
other situations, e.g., an open borehole or a cased borehole. The near rock earth 
formation 305 may include a gravel pack. 

[0028] Referring now to Fig, 5, the nodal configuration for the model of Fig. 4 is 
shown. A different color is used for each of the four regions. After application of the 
boundary loads and/or displacements the result is the nodal configuration shown in 
Fig. 6. An alternate image that may be obtained is the one shown in Fig. 7 wherein 
the bonds that have been fractured during the deformation process are shown in black. 
The failure of the borehole is easier to interpret than in the illustration of Fig. 6. 

[0029] Another capability of the present invention is demonstrated in Fig. 8. The 
illustration on the right corresponds to a fracture simulation wherein the near rock 
compressibility is 1 .2 times the compressibility of the far rock whereas the figure on 
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the right corresponds to a situation in which the near rock compressibility is 0,8 times 
the compressibility of the far rock. 

[0030] Those versed in the art would recognize that in a vertical well, the stresses and 
5 deformation are generally azimuthally symmetric whereas in a deviated well, the 

stresses and deformation will be azimuthally asymmetric due to the radial component 
of the overburden stresses. In the present invention, this possible azimuthal 
asymmetry is taken into consideration. 

10 [0031] While the foregoing disclosure is directed to the preferred embodiments of the 

invention, various modifications will be apparent to those skilled in the art. It is 
intended that all variations within the scope and spirit of the appended claims be 
embraced by the foregoing disclosure. 
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